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Abstract: We consider a finite quantum system coupled to quasifree thermal reservoirs at different temperatures. 
^ Under the assumptions of small coupling and exponential decay of the reservoir correlation function, the large 
deviation generating function of energy transport into the reservoirs is shown to be analytic on a bounded set. Our 
method is different from the spectral deformation technique which was employed recently in the study of spin- 
boson-like models. As a corollary, we derive the Gallavotti-Cohen fluctuation relation for the entropy production 
and a central limit theorem for energy transport. 
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1 Introduction 



1.1 Fluctuations in open quantum systems 

^ \ 

^ Recently, the physics community has shown quite some interest in current fluctuations in nonequi- 

Tj- ■ librium quantum systems. We mention two interesting perspectives: 

■ 

' 1) Since the work of 111411171 , it has become clear that nonequilibrium systems, both classical 

■ and quantum, exhibit a symmetry in the fluctuations of entropy production. This symme- 
^ . try, dubbed the "Gallavotti-Cohen Fluctuation Theorem" holds far for equilibrium. 

^1 2) It has been realized ||29l that noise between electron contacts shows distinct signs of Fermi 

■ statistics, studies of this kind go by the name of "Full counting statistics". 

?H \ Perhaps the most important promise of fluctuation theory is in the construction of nonequi- 

librium statistical mechanics: Via the study of the large deviation rate function, one hopes to 
find a useful variational principle describing nonequilibrium stationary states. Recent papers 
taking part in this project are e.g. l|4l l3n[T2| . 

In this paper, we study heat current fluctuations in a nonequilibrium model of the type 
'spin-boson'. We prove that the large deviation generating function corresponding to energy 
transport exists in a bounded (but arbitrarily large) set around and that it is analytic. 



"^Postdoctoral Fellow FWO-Flanders, email: wojciech . deroeck@ f ys . kuleuven . be 
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1.2 Large deviation generating function 

We briefly sketch the framework of large deviations. 

Assume that we have a family of M'^-valued random variables At, indexed by time t E 
and with distribution given by the expectation E^. To fix thoughts, one can think of the At as 
time-integrals of some variable a(s), < s < t, i.e. 

Af = [ dsa{s) (1.1) 
Jo 

The large deviation generating function on (if it exists) is defined as 

/(k) :=lim^logEt(e-('^l^*)) (1.2) 

where (■ |-) is the canonical scalar product on M.^. 

From the function /(k), one can extract large time properties of the observables At, as sketched 
below, see [6J for precise statements and details. 

1) If / is analytic in a neighbourhood of 0, then At satisfies a central limit theorem with mean 
— 1^(0) and covariance a = |^(0) (the gradient and the Hessian of /). Let 



then 



limEt(e-^(^l^')) =e-(^l"^), 7 G M'^ (1.4) 



2) If / is differentiable on M.^, then the family ^ satisfies a large deviation principle with rate 
function 1(a) given by 

/(a) :=- inf + /(/?)) (1.5) 

Heuristically, this means that 

Probt(y ^a)~e-*^("\ aeR'^,t^oo (1.6) 

(in a logarithmic sense) as t t 00. 

In classical statistical mechanics, the existence of the large deviation generating function can 
usually be established through a convexity argument, see e.g. ||37| . A similar general under- 
standing is lacking in quantum statistical mechanics (see however [1} [19l |28l |33l for partial re- 
sults). Another -even conceptual- problem in quantum statistical mechanics, is how to describe 
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joint large deviations of several noncommuting variables. Remark that it was exactly to solve 
such a conceptual problem for the central limit theorem, that the framework of the fluctuation 
algebra was constructed flSj. 

We consider a quantum setup where At corresponds to the total heat transport into reser- 
voirs. Hence the setup is somewhat different from that in pTKTlUlSHlSl; the expectation E{g{At)) 
for some function g can not be formulated as an expectation of some observable in a quantum 
state, rather it is the probability of obtaining certain (differences of) measurement outcomes. 
The problem of joint distributions for non-commuting observables does not even appear in this 
context since the different reservoir Hamiltonians do mutually commute. This is discussed more 
extensively in [11 J. Our result will establish the existence and analyticity of /(k) on a compact 
(but arbitrarily large) set containing 0. Hence, we do not prove the large deviation principle (but 
we do prove the central limit theorem). The Gallavotti-Cohen fluctuation theorem is a simple 
corollary of our result. 

1.3 Open quantum systems with finite reservoirs 

Our model describes a small quantum system (an atom) interacting with a quantum system 
with many degrees of freedom (a reservoir). We choose the reservoir as simple as possible: a 
free field of bosons, although fermions would do just as weljE The system is coupled to the 
reservoirs through a term, which is linear in the field creation and annihilation operators. This 
type of models are known as Pauli-Fierz models, or, in the simplest case, the spin-boson model. 
These models arise as toy-models in solid state physics, were the bosons are lattice phonons, 
or through the dipole approximation in QED, where the bosons are photons, see jZl for more 
background. 

To make the statements mathematically sharp, we consider this field in the thermodynamic 
limit, or equivalently, in the limit where the modes form a continuum. However, for the sake 
of distilling the right physical question addressed in this paper, we start from a finite-volume 
setup. 

1.3.1 Setup 

Fix a finite-dimensional Hilbert space £ with self-adjoint Hamiltonian E and let /C be a finite 
set which indexes the heat reservoirs at inverse temperatures PkeJC > 0. The superscript n G N 
indicates that the thermodynamic limit {n / oo) has not yet been taken. See also Section [T74l for 
specific notation and conventions. To each k E K,, we associate 

1) A finite-dimensional one-particle Hilbert space ^k,n and its bosonical second quantization 

^s{^k,n)- 

^In fact, they would simplify the technical work 
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2) The coupling operator V^.n G B{S, E ® f)fc,n)- 

3) A self-adjoint one-particle Hamiltonian hk,n acting on with corresponding second 
quantization dr(/iA; „). 

4) A Gibbs state Pfc,/3fe,n on i3(rs(t)fc,n)) at inverse temperature 

PK,.AR\ - ^^^,-,..r(...„)]^ fleB(r.(^u) (1.7) 

We define the total interacting Hamiltonian on E ®k&K ^s{^k,n) as 

Hx,n = + J]dr(/i,,„) + Aj](a*(yfe,„) + a(yfc,„)). (1.8) 
keK. keic 

We take as initial state 

Pf®PR„, Pr„ := ® PfcA.n (1-9) 

corresponding to initially decorrelated reservoirs and an arbitary state ps on B{E). 



1.3.2 Transport fluctuations and their limits 

We introduced the finite volume systems in order to pick the right expression for transport 
fluctuations, and hence, now that all tools are in place, we ask what we mean by transport 
fluctuations in the finite-volume models. 

Note that the reservoir Hamiltonians dr{hk,n) mutually commute and that they have dis- 
crete spectrum. Hence one can measure them simultaneously in the beginning and at the end of 
an experiment. To determine the transport (of energy), we look at the differences of those mea- 
surement values. Let T := Hfee/c sp(dr(/ife „)) and let P^er be the joint spectral projections of 
dr(/ijfc „) corresponding to the eigenvalues x — {xk)keic- The standard interpretation of quantum 
mechanics yields the probabilities 

^Pe,t,xAy) E Ps®PKn[Pxe-''''-''P.'e''''"'^P.] (1.10) 

X ,x' £T ,x' —x=y 

for observing energy differences y G R '^L The Fourier-Laplace transform of this measure has 
a nice expression which is better suited for taking the thermodynamic limit: Using that (the 
density matrix corresponding to) commutes with the spectral projections P^, one arrives at 

/ dy P,„t,A,n(y)e-('^l^) ^P£0 PR„ [r(^_,,„)e'*^^-r(^/;,,„)e-*^^-] (1.11) 
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where 



w.,n = ( © e-'^''^'^-) (1.12) 

keK 



We will study the infinite-volume limit of this expression, given by (|2.9|) and introduced in 



Section I2.2[ In Section |3] we will substantiate the claim that (|2.9|) is indeed the n | oo-limit of 



dLTTj) 



This approach to fluctuations was already used in [|27l |34l l24l |38B for fluctuations of heat 
and work, and, most widespread, in [29il30| for fluctuations of charge transport ("Full counting 
statistics"), made mathematically transparant in 



1.4 Conventions and Notation 

For £ a Hilbert space, we use the standard notation for 1 < p < oo 

BpiS) := {S e B{8), Tr [{S*SY'^] < oo} (1.13) 

and 

:= (Tr (1.14) 

For a Hilbert space f) we write 

r,"(f)) := Sym, ®" f), V,{^) := © r:(f3) (1.15) 

neN 

where Sym„ projects on the fully symmetrized subspace and Ts{^) is the bosonic Fock space 
built on f). For operators C on f), we write (whenever the RHS is well-defined as an operator on 

r(c) = 0®"c (1.16) 

n 

c^dr(c) = 0Vi®...i® c ®i...®i (1.17) 

„gN i=l i'th position 

For W G B{8,£ ® f}), we use the generalized creation and annihilation operators a{W) / a{W*) 
on£ ® rs(f)), (see IZl for an extensive review of this notation). If, for some G [) and D e B{£), 

Wu = Du0tp, ue£ (1.18) 

then a*{W) = Cg> where a*(^/') is the more familiar creation operator . 

For a Hilbert space P), we write f) for its conjugate space, which is fixed by an antiunitary 
map [)^P):a(-^a. If[) = L'^{X, C) for some measure space X, the map a i-^ a is identified 
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with the complex conjugation on functions A* — > C. If i? G then R G B{1)) is defined by 

Ra = Ra. 

For K E C^, we write 

= (JR,Ki, . . . , ^Kd), = i^f^l, ■ ■ ■ , S/^d) (1.19) 

For indicator functions, we use the notation Ind(-), i.e. for a premise a{x) dependent on some 
variable x 

J Ar f \\ / 1 a(x)istrue , , 

Ind(a(x)) = I Q «(^) is false ^^'^^^ 

1.5 Outline 

We introduce the model in abstract terms in Sections 12.11 and 1 2.21 immediately followed by the 
result in Section |Z3| The physical justification of this model is given in Section |3^ where it is 
explained how it emerges from the quantities discussed in Sections 11.3.11 and ll.3.2[ In Section 
13.31 we discuss related results in the literature. The rest of the paper is devoted to the proofs. 
Section IHcontains the main line of reasoning and the technical lemma's are postponed to Section 
|5l The main idea is in Lemma l4!4l whose main ingredient is Lemma |4!2| 



2 Model and results 
2.1 Zero-temperature objects 

Introduce a finite-dimensional space £ with a self-adjoint Hamiltonian E and (for each k e 
/C) one-particle spaces f)/,. oo with a self-adjoint operator /i^ oo on i)k,oo- We also need coupling 
operators Vk^oo e B{£,£ f)fc,oo)- 

One should think of these objects as defining the zero-temperature Hamiltonian of the sub- 
system+reservoir system, formally 

Hx,o. ■.= E + J2 dT{hk,oo) + XY1 («(^fc,°o) + a*(\4,oo)) (2.1) 

The heavy notation with the subscript oo is because in what follows, more natural infinite- 
volume objects are introduced. The objects with subscript oo are relevant at /3 = oo. 
We anticipate the finite-temperature by introducing the Bose-density operators 

a,oo := (6^'="'=- - 1)-^ (2.2) 

We will use the above notation to build a dynamical system which represents our system at 
positive temperature. The connection between the finite-volume objects, introduced in Section 
|1.3[ and the inifinite-volume model, is given in Section |3] 
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2.2 Positive temperatures 

Define 



f)fc := t)fc,oo © [)fc,oo f) := © f)fe 
fce/c 

^fc := hk^oo © (-/ifc.oo) := © /ifc 

^fc := v^l + a,ooVfc,oo © ^/C^Vk,oo V := ® Vk 
Wk,Kk '■= e"'"='' := © WfcAfc 

fcg/C 

Let the total Hilbert space be H := £0T^{[)) and define on £ (g) {V{dT{h))nV{a{V) + a* (V))) , 

Hx:=E + dT{h) + X{a{V) + a*{V)) (2.4) 



The following theorem comes from IHOII 



Theorem 2.1. Assume that \\V\\ < oo. Let Hq := E + dT{h) and denote 

I[u) ■= e™^«(a(V) + a*(\/))e-'"^°. (2.5) 

The series 

jjX.^^itHoY^ j dMi...dM„ /(«„)... /(Ml), (2.6) 

"'^^ 0<Ml. ..«„<< 

originally defined on the dense subspace Vi (see Section \4.1h extends to a strongly continuous unitary 
group on H, which will also be denoted U^. Its self-adjoint generator is an extension of Hx as in (|2.4|) 
and will be simply called Hx in what follows. 

Let p£ be a state on B{£) and let pr be the state on B{Ts{1:))) given by 

PK[-] = {n,-n) (2.7) 

where Vl = 1©0©0... G rs(f))is the vacuum vector. We will take p£ © pr as initial state on 
B{H) for our dynamics. Unless otherwise stated, we assume ps to be arbitrary. 
We now introduce our main object of study 

Assumption A-1 (Bounded interaction). For all k with 3?/t G D, 

\\wiiV\\ < oo, < oc (2.8) 

The following lemma follows from Section |4H 
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Lemma 2.2. Assume there is an open set D c M''^' with G -D such that Assumption (|A-1|) is satisfied 
and let be as in Theorem IZll Than the function 

K^pe® PR [T{w-,)U^j:{w^)U^] (2.9) 

/ifls an analytical continuation from {^k = 0} into {^k g -D}. 

The function (|2.9|) should be thought of as the Fourier-Laplace transform of the probability 
distribution of energy transport. This is discussed and justified in Section |3l 

2.3 Results 

To continue, we need additional assumptions. The next assumption basically establishes that 
the operator hon[) has absolutely continuous spectrum. 

Assumption A-2. There are measure spaces {X, dx) and {y, dy) such that 1} = L'^{X, dx), X = y 
and dx = dyd^ where d^ is the Lesbegue measure on R. For {y, ^) = x e X,we write ^{x) = ^for the 
projection on M. The operator h acts by multiplication with ^{x), 

{h-^){x) = ^{x)iIj{x), e f) (2.10) 

Remark that one can associate to a measurable function X 13{£), which we denote 
X i-^ V{x) and which satisfies 



{v iS)'4),Vu)£^t)= / dx^p{x){v,V{x)u)£, u,v e £,ip e \:) (2.11) 

Define the reservoir time-correlation function 

p.{t):= sup \\V*Se-'"'w,V\\ (2.12) 

Sel3{£),\\S\\=l 

Assumption A-3 (Decay of bath correlations). Let be as defined above. There are C,a > such 
that for K e D 

< Ce-°l*l, tGM (2.13) 

Introduce the set of Bohr frequencies JF := spE — spE and let Is^ stand for the spectral 
projection of on e G spE. 

The following assumption expresses that the coupling between system and reservoir is suf- 
ficiently effective. 

Assumption A-4 (Fermi Golden Rule). 
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1. For all uj E T and dy-almost all y E y, the function V : X ^ ^{^) i-S continuous on the set 
{x = (y, 01 ^ = T'^^'s i^lpliss that V{x = {y, uo)) is well-defined. 



2. IfSe B{£) satisfies 



E E j '^y\\ uy{y,uj)uA \\ = o (2.i4) 



oj — e — e 



then S = clfor some c G C. 

Now comes our main theorem 

Theorem 2.3. Assume Assumptions \A-1\ \A-2[ \A-3\ and \A-4\ There is a Xq > such that for X e 
[— Ao, Aq] and k E D, 

A) := \imt-Hogpe®PK [T{w^^)U^_,T{w^)U^] (2.15) 

exists, is independent of ps and real-analytic in k and X. 

As the main corrolary, we state the Gallavotti-Cohen fluctuation theorem for the entropy 
production. This requires an additional assumption 

Assumption A-5 (Time-reversal invariance). There is a anti-unitary & onH such that for all A G M, 

e-^Hxe = Hx, e~Mr(/ifc)e = dr(/ifc), e{s ® n) = {s n) (2.I6) 

and 6 is an involution, i.e. = 6. 

Theorem 2.4. Assume \A-1\ \A-2\ \A-3\ \A-4\ and \A-5\ Let f{K, A) be as in Theorem |2J] and define for 

= z/(/3i,/32,...,/3|k;|)gMI^' (2.17) 
For z/ G M such that k{1 — u) G D, 

fiKiu),X) = fiKil-u),X) (2.18) 

By Bochner's theorem, there is a nonnegative Borel measure clPp^ < a on RI'^I such that 

Pe ® PR [riw^^)U^,T{w^)U,'] = [ dP,,,i.A(y)e-(''l^) (2.19) 
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for 3ftK G D. Putting k = 0, one sees that dPp^^t^A is a probability measure. It is the infinite- 
volume analogue of the probabilities Pp£,i,A,n introduced in Section [l .3.21 We write Ep^ t^A[-] for 
the expectation w.r.t. (Wp^^t,\- 

The M'^-valued random variable y = {yk) is interpreted as the energy transport into the dis- 
tinct reservoirs. Remark that in thermodynamics, one interpretes S := Y^keic f^^Uk as the entropy 
production. 

Since 

/(«:(//), A) = Wm- log^p^^tAe-"'] (2.20) 

one sees that f{K,{iy), A) is indeed related to (large) fluctuations of the entropy production. 
The following corollary follows from Theorem 12. 3l by [5J. 

Corollary 2.5. Assume the assumptions ofTheorem \2.3\ Then the W^-valued random variable y satisfies 
a central limit theorem with mean — ^/(/t, A) and covariance ax ■= ^/(k, A) Let 

then 

Ep,,t,A[e-'(^'^*^] — .e-(^l"^^\ 7GMI'^I (2.22) 

The expectation value — ^/(k, A)|^^q and the covariance ^/(k, A)|^^q can be written in a 
more familiar form. Introduce the operators 



Ak,t := UhA^ihk)U,^ - dVikk), (2.23) 



Then 



-^/(/^,A)L^o = limip^®PR[Afc,J=:(A,) (2.24) 
g^^f{f^,\)\,^, = hm-pe(E)pK[{Ak,t-t{Ak)){Ak>,t-t{Ak>))] (2.25) 

where the convergence of the expressions on the RHS is a consequence of the analyticity of 
/(k, A). However, it is not true in general (beyond second order in k) that 

Pe ® PR [T{w.^)U^,r{w^)U,^] =pe® pR[e-^^'^'=^'=.*]. (2.26) 

See [llj for a thorough discussion of different approaches to quantum fluctuations. 
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3 Discussion 



3.1 Initial state 

We formulate our result only for particular intitial states, namely ps ® Pr with p^ the vacuum 
state. One could ask whether Theorem 12.31 still holds for a different inital state. In fact, by a 
slight generalization of our method, one can prove (see e.g. the previous version of the present 
paper) that the same result holds if one replaces 

P£ ® PR [T{w.^)U^,r{w^)U,^] (3.1) 

by 

P8 ® PR [u^Mw^.)u^tny^.)u,'u^] (3.2) 

for arbitrary s. That is, /(/t, A) is independent of s. 

However, in Section [1.3.21 one sees that the very choice of our object of study (|3.1|) depends 
on the fact that pr is 'diagonal' in the operators dT{hk). This (or rather, its finite-volume ana- 
logue) is used in going from (|1.10[) to (|1.11|) . Expressed more dramatically, an expression like 
(|3.2|) does not appear! 

3.2 Thermodynamic limit 

We skipped over a thorough justification of the object (|2.9|) , which features in our results. We 
remedy this by telling in which sense the dynamical system is the infinite-volume version of the 
finite-volume systems and how the expression (|2.9[) emerges. Usually, thermodynamical limits 
are constructed by specifying volumes which go to infinity in some sense (e.g. in the sense of 
Van Hove). In our case, such an explicit setup is not necessary (though of course possible). 
We simply demand the following relation between the finite-volume objects and the objects 
introduced in Section IZTl 

Assumption A-6 (Thermodynamic limit of finite-volume models). Let 

9iA^) = ^^,x _ ^2,t(x)=e-'*^(l-^^^^), tGM,xeM+ (3.3) 



For S G B{S), i = 1,2, we have 
and 



\\Vl^SgiAhk,oo)Vk,oo\\ < oo (3.4) 

yk,nS9iAKn)yk,n — ^ VT^^^t (/lfc,oo) V'fc.oo (3.5) 



n oo 



uniformly on compacts in t G 
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If Assumption I A-6l is satisfied, a large class of correlation functions converges. There is quite 
some arbitrariness in this statement, which is usually not considered in the literature. 
Define 

^kAt) ■■= e'*^°-(a(\4,„) + a*(rfc,„))e-^*^»- (3.6) 
<l>fc(t) ;= e'*^°(a(\4) + a*(Vl,))e-'*^° (3.7) 

(3.8) 

Assume Assumption |A-61 then for alH, t' e M and S eB{£), 

p£®UJ^,n[^k,n{t)S^k',n{t')\ p e ® Pk[^ kif) S ^ y it' )] (3.9) 

Of course, from p.9|) one deduces also convergence of higher-order correlation functions (since 
the states c<jr „ and pr are quasifree, those are expressed in terms of the second order correlation 
function). In particular, one has also convergence of the same correlation functions with the 
time-dependence now given by the fully interacting evolution, that is, let 

:= e^*^^^-"e-"^°'"$fc,„(t)e^*^«'"e-^*^^''' (3.10) 

;= e"^^e-"^«$fc(t)e'*^°e-"^\ (3.11) 

then equation (|3.9|) holds with $^ replacing $, as follows from a Dyson expansion, e.g. (|4.3|) . 
It is now straightforward to see that Assumptions lA-l] and lA^ imply 

pe ® Pr„ [r(w_,,„)e'*^->'T(z/;,,„)e-^*^^-] ps ® Pk [T{w^^)U^_J:{w^)U^] (3.12) 

ntoo 

where the LHS was introduced through physical considerations in Section [l.3.2[ 

The critical reader might wonder why there is in our presentation no mention of W^*-algebra's, 
which often play a promiment role in the mathematical formulation of statistical mechanics. If 
one defines the Araki- Woods algebra 21 as in Section 14.61 one finds that the dynamics 

2t 9 A U\AU^, s G M (3.13) 

leaves 21 invariant. Physically, one should restrict the state pc ® Pr, originally defined on B{£ ® 
rs(P))), to 2t. However, in our approach, it is neither mathematically nor physically necessary 
to consider this restriction. We study the expression (|3.1|) , which is well-defined and whose 
motivation is via (|3.12|) . 

For the same reasons, we do not have to ask ourselves whether the operator (|2.4|) is the right 
choice. In the literature, this operator is called the semi-standard Liouvillean, but one can also 
consider the standard Liouvillean. Again, the resolution of any possible ambiguity is via finite- 
volume limits. That being said, it might be worth remarking that (|3.1|) can be expressed as the 
expectation of powers of a relative modular operator, see H32l , thus providing a more algebraic 
starting point for our work. Another possible approach is in [2], where the expression (|3.1|) is 
constructed (for fermions) via different, but essentially equivalent reasoning. 
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3.3 Comparison with other works 

There has lately been a lot of work on spin-boson and spin-fermion models, or more general, 
Pauli-Fierz models. 

We feel our work is technically closest to |22|, in which one considers the spin-boson model 
and one proves that the generator of the dynamics has absolutely continuous spectrum for A 7^ 
0, except for one eigenvalue which corresponds to the stationary state. The other eigenvalues of 
the system at A = turn into resonances whose location is in first nonvanishing order predicted 
by the Lindblad generator. The assumptions are very similar; to allow for a comparison, we 
assume that 

Vu = Du^ip (3.14) 

for some D e B{S) and ?/> G [) ~ L'^{R, L'^{y, dy)), in which case a*{V) =D®a*{'4)). 
The basic assumption in ||22]| reads 

Assumption A-7 (analytic coupling). The function i/j is analytic in a strip {^z < 6} and 



sup / de||^(e + i7)r <oo (3.15) 



Assumption I A-7l implies that 



< Ce''^ (3.16) 



which is just Assumption |A-3| for k = 0. However, the k have no analogue in f22| and we would 
need to assume Assumption I A-7| with ip derived through (|3.14|) from wiiV rather than from V. 

In contrast, we do not need any additional infrared condition on ^ i-^ ^"(0/ contrary to 1122 II . 
This is because we construct the dynamics via the Dyson expansion instead of via the Nelson 
commutator theorem. Physically speaking there is of course already an infrared condition 
present since 

ll^ll <oo ^ Yl \\iPkhk,oo)-'^'Vk,oo\\ < 00 (3.17) 
fce/c 

with the notation as in Section IZTI 

The technique of ||22|| consists of a spectral deformation of the generator Hx. We employ 
time-dependent perturbation theory and we rewrite the Dyson expansion as a one-dimensional 
polymer model. This is embodied in Lemma |42l Starting from that lemma, one can obtain our 
result through a simple cluster expansion (as in the previous version of this paper). However, 
since the polymer model is one-dimensional, we can apply the transfer-matrix technique. In 



^That is, in terms of the zero-temperature coupling operator, or 'form-factor' T4 , 
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dealing with the transfer matrix, we use a variant of the spectral deformation technique, such 
that our technique is not as different from [|22|| as might seem. 

Assumption IA-31 cannot be weakened without changing the method drastically. Note that 
one cannot assume Assumption IA-31 for D = M'^I since that would imply that 

rI^I 9 K h^p^(t) (3.18) 

is a bounded analytic function, hence constant. 

Results that need weaker regularity properties of V^(^) are e.g. BH, EUHl, flSl. In those works 
one employs Mourre theory or renormalization group techniques, however they do not permit 
to localize the resonances. 

A different type of works are those using scattering theory. This approach was initiated in 
Il36ll , but so far, it has not been successful for spin-boson type models, although it works well for 
junctions l|T6ll . 

From the physical point of view, our result is closer to ||23l |20B where one studies a non- 
equilibrium setup and one derives approach to a non-equilibrium steady state and the Green- 
Kubo relations, or to [[2T|, where one studies a form of the central limit theorem. See |[TT| for 
an extensive discussion of the difference and similarities of different approaches to quantum 
fluctuations and central limits. 



4 Proof of Theorems 2.3 and 2.4 



4.1 Construction of the dynamics 



Let In be the projector on S ® Tl'{i)) (the n-particle sector, see Section n~4|) and let the domain 

Vi c £ ®T,{[))he defined by 



(4.1) 



Letifo := E + dT{h) and 



(4.2) 



For = 0, the series 



E 



dui . . .du„/«(u„) . . 



(4.3) 



0<Ml...M„<t 
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(4.4) 



Since the argument in [jTO| showing that (|4.3|) is a a strongly continuous group on Vi, depends 
only on the assumption < oo, this remains true for k satisfying Assumption lA-ll and (|4.3)) 
can be taken as the definition of (|4.4|) . 

In what follows and unless stated otherwise, we will assume that Assumptions lA-ll and lA^ 
are satisfied and that G D. 

4.2 Dynamics and notation on ;Bi (7-^) 

It is advantageous to rewrite the object of study in a slightly more abstract way Let Pi,® stand 
for the subspace oi Bi{£ ^TgCi))) defined by finite linear combinations of |0i)(02| for0i,02 £ T^i- 
From the conclusions of Section 14. 1[ it follows that 



)) =: Zt\A) 



(4.5) 



maps T>i^^ into itself. In what follows, we write 



M{S) := i[E,S] 



(4.6) 



as a bounded operator on B{£). We define the embedding T| : B{£) ^ Bi{£ ^s{^)) 



(4.7) 



and the compression Jj : Bi{£ ® rs({))) — > B{£) 



(4.8) 



with Trr^((,) the trace on Bi(TTr^{t))) (Hence is actually a partial trace). 
We have hence rewritten 



(4.9) 



where ps is the density matrix, corresponding to the state pe, i.e. ps [S] = Tr [psS]. 
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4.3 The deformed Lindblad generator 

If \\Pk\\i '■= J dtpK(t) < oo for K = 0, we can define 

POO 

T = ~iJ2 E / Uy*le,e-''^''-^^Vl£At. (4.10) 



e—e—Lxj 



Assuming additionally the first statement of Assumption I A-4[ we introduce the deformed Lind- 
blad generator. For S G let 

L^{S) = -i{TS- ST*) + 271 J2 E ^£.V*le^,S5{h-u;)w^Vle^ (4.11) 

ojS^ e—e'=LU 

where the operator-valued Dirac-delta distribution 5{-) is well-defined by the continuity as- 
sumption in Assumption I A-4l For example, one can take a sequence of functions converging in 
the sense of distributions to 6{- — u), then the mentioned continuity assumption assures conver- 
gence in (|4.11[) One checks, see e.g. fTO\, that for k = 0, or equivalently, = 1, we recover the 
usual definition for the Lindblad generator, which satisfies 

Tr[L,=o(^)] = (4.12) 

However, since the second term in (|4.11|) is a completely positive map, it follows that e*^" is a 
completely positive semigroup for {3?k G D}. 
We need the following properties of L^. 

Theorem 4.1. Let be as in (|4.11|) and M as defined in Section \4.2\ 

1) Assume Assumption \A-4\ and fix a t < 0. The operator e'^^" has a maximal simple eigenvalue e^^" 
with /(k) g M and there is a 'gap' > such that 

sup{|;z|,;z G sp (e"^"^ ) \ e"-^'= } < e"-^'=(l - e-"^'^) (4.13) 

The eigenvector corresponding to e"^-^" can he chosen a positive invertible operator. 

2) 

[L.,M]=0 (4.14) 

3) Assume < oo. For all t > 0, 

II J|Z^iXXt - e-'"(^''^^+'^«) II ^ (4.15) 

where the LHS is continuous in A, k, t. 

Statement (1) of Theorem 14.11 is the only place where we use the second statement of As- 
sumption |A-4[ It is a non-degeneracy assumption which enters the non-commutative Perron- 
Frobenius theorem. 
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4.4 Dyson expansion and transfer operator 

Our basic tool is a rearranged Dyson expansion, whose properties are collected in the upcoming 
Lemma |4!2| Fix a parameter r > and define onB{£) for n E No, 

W^'^'- := Z^,X (1 - W •••(!- m Zt'X (1 - W Z^'X (4-16) 

n-1 2 1 



(n — 1 factors of (1 — X|X|) inserted). The definition (|4.16|) makes sense since Zf'^ maps Pi ,g, into 
itself (see Section 141]) and, obviously, X^B{£) C Vi^^. Whenever reasonable, we will abbreviate 

Lemma 4.2. Lei >V„ = W^'^'^ he as above. 

1) for all m G Nq, 

XiZ,^d_XT = 5^ ^ W„,...W„,W., (4.17) 

2) Assume Assumption \A-3\ There is c := c(/t, A, r) > 0, vanishing as X [ and continuous in the 
three parameters, such that for n> 1, 

llWnll < c^-' (4.18) 



In what follows, we use the Hilbert space /2(No) ® B2{£). Let for n e No, e„ be the canonical 
n'th base vector in /2(No) and let S be the unilateral shift, defined by (setting cq := 0) 

Sen = en-i (4.19) 

Recall that £ is finite-dimensional, which allows to define the embedding P„ : B{£) /2(No) ® 

B2{£) : u ^ Cn® u and compression P* : Cn ® u ^ u. We are led to examine the following 
operator on /2(No) ® B2{£); 

rpn,X,r ^ P„WnPl +S®1 (4.20) 

neNo 

From Lemma l42l l), one has 

X|Z^;'"^"X| = p*(r'^'^'")'"Pi (4.21) 

If the operator T := T'^'^''^ had a maximal eigenvalue, isolated from the rest of the spectrum, 
we could easily estimate the n y oo asymptotics of (|4.21[) . However, upon realizing that 



sp5 = {zeC, \z\ < 1} (4.22) 

this surely fails at /t = 0, since the highest eigenvalue of e''^'==° is 1. This difficulty is addressed 
in the next section. 
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4.5 Spectral deformation 

Introduce the unbounded operator 

R=Y,nPnP:. (4.23) 

nSNo 

The following statements are straightforward. 

Lemma 4.3. For6eR and W G B{B{£)), 

1) e-^^^e-"^^ = e-^^ 

2) e^^PnWP^e-^^ = e('^-"^)^P„l^P^ 

Most importantly, the operator Q^^Te~^^ does have an isolated eigenvalue for well-chosen S, 
as we show now. 

Lemma 4.4. Let 6 := — 1/2 lnc(K, A, r), where the latter was introduced in Lemma W2\ There is a 
Ao > such that for A G [— Aq, Aq], k e D and r varying in some compact set D^, the operator 

QSRrp^-5R (424) 
has a maximal simple eigenvalue e'^-^"-^-^ with f^^x,T ^ I^- There is gt,xr > such that 

snp{\z\,ze sp(e^^Te-^^) \e"^'=-^--} < e"-^«.^.-(l - e-'^"-^--). (4.25) 

The eigenvector Gk,\,t corresponding to this eigenvalue can be chosen such that P^Gk.a.t G B{S) is an 
invertible, positive operator. The function fK,\,T is real-analytic in k ^ D,\\\< \q and r G D^. 

Proof. By Lemma l43l 
where 

AT := e--^^ + (Wi - Q-^-i^-'M+iL^)-^ ^ J2 e^''-^^^ PnWnPl (4.27) 

n>l 

By Lemma |42l (let cbe as defined therein) and assuming \ce^\ < 1, 

II AT|| < e-^ + ||Wi - e-'-(^"'^+'^'=)|| + (||P*|| sup ||P„||)^^ (4.28) 

neNo 1 - ce^ 
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The norms 1 1 -P* 1 1 , 1 1 Pn \ \ are independent of n and finite since dim E is finite, and hence, using 
Theorem 14. ir 3), || AT|| vanishes as A | and as 5 = 5. 

Remark that M is self-adjoint on B2{S) and that spM = JF. By Theorem I4.ir 2), we can hence 
decompose = (Bc^eJ^L^^^ where L^^o; acts on the cj-eigenspace of M. Hence 

{z - e-^r{>^-''M+iL^)-^ = © e'"^"'" (e'"^"'"z - e"^''^-) (4.29) 



Theorem I4.1l fl), the expression (|4.29|) and compactness of the unit circle in C yield that there is 



a e > 0, C > 0, such that for e"^-^"* — e < \z\ < e"^^", for k E D and for r varying over some compact 
set. 

The existence of an isolated eigenvalue and positivity of the eigenvector now follows from (|4.28|) 
by standard perturbation theory, see e.g. [25J. Positivity of the eigenvalue follows since by (|4.5|) , 
Zf''^ is a completely positive map for = 0. 

Real Analyticity in k and A for A 7^ follows from analyticity of and AT, both of which are 
straightforward consequences of Assumption lA-ll Since e'*'^ doesnot have a limit as A | 0, 
analyticity at A = is not immediate. However, since, /k,a,t is analytic for A 7^ and continuous 
at A = 0, it is analytic. □ 

By Lemma l434l we get for m large enough 

logXi^r/™^^! = Ux,r + log + 0(e--^-^-)) (4.31) 

lilt I 1 1 L 

where Pq^ ^^is the projection on Gk^a.t- 

Taking r, r' G D^. such that mr = m'r' for some m, m' E N, we get from (|4.31|) that ffc,x,T = 
/k^a,t'- Since /K,A,r is also continuous in r, it is constant and we write /^^a := /k,a,t- Theorem 12.31 
now follows with /(k, A) = A^/^^a by (|4!9|) . 



4.6 Proof of Theorem Q 

Assume that 

Let il be the Vr*-algebra (Von Neumann- algebra) which is generated by the sets 

B{S)(^1 and {e^("W+'^*('^»,V^ e f)} (4.33) 
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Remark that the expansion (|4.3|) shows that for all t, 

^tHo^-itH^ e il. (4.34) 

(See HZl for details on iy*-algebra's). Extend the notation k(z/) in Theorem 12.41 to u E C. The 
maps of automorphisms 

il 3 A ^ r]siA) := r(w;,(i,))Ar(w;,(_i,)), s G M (4.35) 

is a Ty *-dynamics and ps ® pr is a 1-KMS state wrt. to this dynamics. This can be easily checked 
or read in the literature, see again |l7|]. Then, the KMS-condition reads that for A, B E il, the 
function 

Pe ® PK[Ar]siB)] (4.36) 

is analytic in {0 < S^s < 1} and satisfies 

Pe ® Pr[Vs{A)B] = pe® p^[Br,,+,{A)] (4.37) 

Choosing A = e~^*^oe'*^^ and B = e~'*^^e^*^°, inserting 1 = 66, using Assumption IA-5[ the 
general property p{C*) = p{C) (true for every state p), [e~'*^°, T{wk)] = and the invariance of 
Ps ® Pr under the dynamics e"^*^" ■ e'*^°, one gets the relation 

p£ ® PR [v~iu{U-t)Ut] =Pe®pK [ri-i{i-u){U-t)Ut] (4.38) 

for — 1 < z/ < 0. This is extended by analyticity to values of u such that k(z/) g D. Theorem 12.41 
follows since by Theorem |2.3[ /(/«, A) is indepenent of p£. 



5 Proof of some estimates 

We prove the lemma's that were used in Section HI As in Section HI we always assume As- 
sumptions |A3] and |A21 and we take k such that {^k E D} where D is as in Assumption lA-ll 



5.1 The Wick-representation of the dynamics on B{S) 

The aim of this section is to introduce a convenient notation to handle the Wick-ordered Dyson 
expansion, stated in (|5.5ti5.6|) . The result is equation (|5.12|) . 

Recall the representation of as a function V : X ^ ^{^)r introduced in l2.11[ Denote 

V*{x) := e"^y#(x)e-"^ teR,xeX, V*{x) = V{x), {V{x)y (5.1) 
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By Assumption IA-2[ both h and Wk can be represented as multiplication operators with func- 
tions on X. We will denote these functions by respectively ^{x) and Wk{x) (consistent with the 
use of ^ in Assumption I A-2[) . 



Introduce the space Z = A" x {1, 2, 3, 4} with elements z = (x, j) and measure dz = dxdj {dj 
stands for the counting measure on {1,2,3,4}) and the maps Q^eK zez ^ ^i^i^))' 



Qu,z={x,j)i^) — ^ 



w^{x) Vu{x) S j = l 

e'"«(^')u7_J(x) V*{x) S J = 2 

e-'"«(^)ti;_J(a;) S j = 3 

ei<{^)wV{x) S V*{x) j = A 



(5.2) 



We now introduce the pairing coefficient C{z, z') for z, z' G Z; 

C{z = x,j- z' = x',j') := 6ix -x')h ^' = ^' { / = 4 ^' = ^' { / = 3 (5.3) 

y otherwise 

For n E 2N, let Pair(ri) denote the set of partitions of {1, . . . , n} in pairs. For such a partition 

71 G Pair(n), we write 

J (i, is one of the pairs in the partition vr , , 

[1,1 ) n ^ < ^/ ^ ^ (5.4) 

The following representation for JiZ^'^J-^ is our starting point. 

ne2N 0<ui<...<u„<t 7rePair(2n) 

A-'^ Jdz,...dzJ n Ciz,, ^^')] Ql-^u„,.„ ■ ■ ■ QUu.,., (5-6) 

It follows from the definition (|4.5|) , the Dyson expansion (|4.3)) and the Wick theorem. 

Let [0, t]2 be the set of (unordered) couples in [0, t] and 

fit := {a C [0,t]2, |a| < oo} (5.7) 

We remark that there is an idenfification between n G 2N, < Mi < . . . < m„ < t, tt G Pair(?2) and 
cr G fit with |cr| = n/2, given by 

a= U {(n^UiO) (5.8) 

(j,j')^7r 
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By writing dn and d„7r for the counting measures on respectively N and Pair(n), we define, 
using the above identification, 

da := dn x dui x . . . x d-u„ x d„7r, (5.9) 
This definition could be ambiguous when |cr| = (hence a = 0), which we fix by defining 

[ dalnd(a = 0) = 1. (5.10) 

Thus, we have made Vlt into a measure space. Using the same identification, we define V'^'^(cr) G 
B{B{£)) to equal the line dS^I 



and we again abbreviate V((t) := V'^'^{a). 
We have hence rewritten (|5.5H5.6|) as 

IlZ^A^^ = e''-''^' [ daV{a) (5.12) 

For convenience, we also define (it C as the set of those a with |cr| = 1. Hence (it is the set of 
ordered pairs in [0, if:]. We will write the elements fo this pair as s((t), with s{a) < s{a). 

We stress that up to this point, nothing happened; we just cooked up a fancy notation, cul- 
minating in equation (|5.12|) , for the Wick-ordered Dyson expansion! 



5.2 Proof of Lemma IQ 

Statement (1) of Lemma |4!2] is an obvious consequence of the definition (|4.16|) , we concentrate 



on Statement (2). We first establish the crude a-priori bound (|5.16|) . 
Let (ua) be a basis in £ and define 



E 



dx{Ua,{V{x)yUa'){Ua",V{x)Ua"')e''^^''^W^{x) 



X 



(5.13) 



Since E is finite-dimensional, the function qj^t) is dominated by a multiple of ^^(t) (as de- 
fined in l2.12|) and vice versa. Using the explicit expression (|5.2[) , (|5.3|) and (|5.11|) , one gets 



) + A ^(gj}K + g_s«:)(- 



S (T - S O" 



=: A~24 



(5.14) 
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One easily checks 

||V(a)||< n X-'d.fl^^^f^) =■■ Gia) (5.15) 

(For example, one can represent V*{x) = Y^aa' Kf ^ (|5.2[) and then factorize 

(|5.11|) ). By a change of integration variables, and summing over all values of \a\, we arrive at 
the a-priori bound 

[ da||V(a)|| <e*ll'='-lli (5.16) 

with II (is 111 = Jf,+ dK{t)dt, which is finite since ||Pk||i is finite. 
Let 

Js,r±{<y) ■■= IndpCT eQt, & C cr, s{a) < s < s{a), s{a) - s{a) ^ r] (5.17) 
One can easily convince oneself that (V stands for the maximum) 



(5.18) 



In words, each a contributing to >V„ contains for each j = 1, . . . ,n — l aa which 'crosses' jr. Or, 
the insertion of 1 — X|Xj forces a pairing to occur. 

Lemma 5.1. Assume Assumption \A-3\ There are c± := c±{k, A, r) vanishing as X I and continuous 
in the three parameters, such that for JL, JTV disjoint subsets of No, 

/da I n J^.r,rMamcr)\\ < (c+)l^+l(c„)l^-l [ daG{a) (5.19) 

Proof Denote by {Gi}i a partition of in subsets Qi satisfying max Qi < min Qi^i for all i. Write 
for the sum over all such partitions. Then, 

LHs of dsn) < 1 5^ n / ( n ^^-.-+ ) ^(^) I (5-20) 

Yl j_ ddJjr,rS^)G{a) \ X daG{a)^ (5.21) 
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Since 



/ daJjr,rA^)Gi^) < ^"^ / dsdsd^{\-\s - s)) (5.22) 

JQt J s<s<s 

< AM duud^{u) =: c_, (5.23) 



hence the first factor in (|5.21|) is bounded by 



\J-\ 



By the argument following (|5.13|) , Assumption IA-31 implies that there are Ck,«k > such 
that d^{t) < C^e-"«l*l (Obviously, ; can be chosen constant if k varies in a bounded set). 
One can bound 



Using ^ . I max Qi — min | > | JV I / we arrive at the upper bound for the RHS of (|5.20|) 



E n / ( n ^^•-.-+(^) ) ^(^) ^ e^l'^+l^e--^-'(V2)|^+|r (^^y^+l (5 25) 



□ 



To conclude the proof of Lemma l42l we use expression (|5.18|l , replacing V — * +, 



< ^ da (^n(^i-,- + Jjr,r+?j < (c- + c+)"-i(e^ll"«ll^)" (5.26) 

To get the last inequality, we represented the product in Ylj=ii'JjT,T^ + JjT,T+) as a sum over 
partitions of{l,...,?2 — l}in2 sets J71 and JV/ we applied Lemma |5J] and we resummed the 
sum over partitions by the binomial formula. Finally, the bound (|5.16|) with t = nr was used. 



5.3 Proof of Theorem D 

These statements are contained in the literature. Statement (1) is a consequence of the Perron- 
Frobenius theorem for completely positive maps, stated in | [T3]1 and valid in our context un- 
der Assumption IA-41 (This is extensively discussed in ||35B ). Statement (2) can be immediately 
checked from the explicit expressions in Section 14.1 1[ For k = 0, Statement (3) is a result of the 
usual weak-coupling theory, see e.g. [lOJ. For k 7^ 0, it is a straightforward generalization of 
these theorems. One can easily follow the arguments in [JOJ and adapt the statements. 
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